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Abstract 

We extend chiral perturbation theory to include linear dependence on the lattice spacing a for 
the Wilson action. The perturbation theory is written as a double expansion in the small quark 
mass m q and lattice spacing a. We present formulae for the mass and decay constant of a flavor- 
non-singlet meson in this scheme to order a and mi. The extension to the partially quenched 
theory is also described. 
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I. INTRODUCTION 



Chiral perturbation theory (x?T) is an important tool for extracting quantitative infor- 
mation from lattice simulations of QCD. The reason for this is that it is impractical to have 
dynamical quarks in simulations that are as light as the up and down quarks, and is 
needed for a controlled, systematic extrapolation in the quark masses. Since describes 
continuum QCD at low-energies, its application in numerical simulations is possible only af- 
ter extrapolating lattice data to the continuum limit where the lattice spacing, a, vanishes. 
In this paper, we study the behavior of the Wilson lattice action close to the continuum by 
incorporating 0(a) effects in a reformulation of ;^PT. A similar approach was first taken in 
Ref. [IJ to investigate the phase diagram for Wilson fermions in two-flavor QCD. 

The quark mass matrix (considering only the 2 or 3 lightest quarks) has a special role 
in QCD — it parameterizes the explicit breaking of the axial symmetries. As a result, the 
light quark masses appear explicitly in the low-energy effective theory. In this paper, we 
exploit the fact that for the Wilson lattice action there is another independent symmetry 
breaking parameter, linear in a ||. To 0(a) this is the only discretization effect, and thus 
a generalization of the chiral Lagrangian can be written which includes all terms linear in a. 

II. EFFECTIVE LAGRANGIAN 

The Wilson action for fermions is given by 

g(W) _ ^ilj(x) / y tl A f j i (x)ilj(x) + ip(x)m q tl)(x) + artp(x)A 2 ilj(x)~\ , (1) 

X 

where 

A^ij){x) = — \U^{x)ip{x + afi) — U^(x — afi)ip{x — a/t)] . (2) 

U^(x) is the gauge group valued field defined on the links of the lattice. In typical lattice 
simulations, r = 1 but we will keep it more general for now. Since we are interested in a 
perturbative study of discretization effects, we follow and consider an effective action in 
the continuum which describes the same physics as the discrete lattice action (including the 
gauge action), well below the cutoff 1/a. The effective action is expanded in powers of a: 

^eff = So + aSi + a 2 S 2 H . (3) 

By construction, Sq is the QCD action. Symmetry considerations restrict the number of 
mass- dimension 5 operators that appear in Si. The equations of motion can be used to 
further reduce the list of operators. One then identifies the operators that already appear 
in 5*0, which give rise to the renormalization of the quark masses and the coupling constant. 
Finally, one is left with a single new term at 0(a), the Pauli term 0: 

aSi = acsw^o-^ V F^ v tj) . (4) 

csw is a constant of 0(l) which is a complicated function of the gauge coupling and r. 

Considering So + aSi as an underlying theory, it has an SU(Nf)L x SU(Nf)ji flavor 
symmetry which is broken down to the vector part by the mass and Pauli terms with 
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coefficients m q and acsw respectively, ipip and ipa^ V F^ v ip break the axial symmetry in the 
same way and therefore, from a spurion analysis point of view, m q and acsw are on equal 
footing, as noted in Ref. M. In particular, acsw is treated as a non-trivial matrix in flavor 
space, just like m q . It is possible to give acsw a flavor structure in simulations by promoting 
the constant r to a matrix. This can supply extra "knobs" which can aid the continuum 
extrapolations. The squares of the Goldstone boson (GB) masses are linear in the symmetry 
breaking matrices of parameters [TJ], conveniently written as: 

X = 2B m q , p = 2W ac S w ■ (5) 

Bq and Wo are unknown dimensionful parameters that appear in the effective Lagrangian 
at leading order (LO) which is defined below. 

%PT for the Wilson action (which will be denoted by W%PT) is an expansion in the 
squares of the small momenta p and the small pseudo Goldstone boson masses. We formally 
consider the expansion to be in two independent small parameters: 

e ~A^~7^ and 5 ~A^ (6) 

X X X 

A x is the scale where new high-energy physics enters and the effective field theory (EFT) 
no longer describes the correct physics. This happens around the mass of the rho meson, or 
roughly A x ~ lGeV. From dimensional analysis, and the fact that B Q and Wq depend only 
on the high-energy details, one can check that the expansion is in fact in m q /A x and in aA x . 

Since Eq. (|3|) is truncated at 0(a) it makes no sense to go beyond 0(8) in W%PT (we 
remark on 0(a 2 ) corrections in the next section). For convenience, we choose to collect 
terms in the LO and next-to-leading order (NLO) Lagrangians as follows: 

LO: £ 2 ~0(e,<5) , (7) 

NLO: £ 4 ~0(e 2 ,e5) • (8) 

The underlying hierarchy consistent with this choice is {e, 6} ^> {e 2 ,e5} ^> 5 2 , and the last 
inequality also implies e 5. In practice, the perturbative expansion should be organized 
according to the actual sizes of the expansion parameters, which are determined by the quark 
masses in the simulation and the size of the lattice spacing. For example, if the simulations 
are done so close to the continuum that S is very small, it might make more sense to have 
the LO Lagrangian be 0(e), and at NLO 0(e 2 ,5). With our convention, 

C 2 = £ tr(0EflEt) - ^ tr(( X + p)£+ + E( X f + p f )) , (9) 

is the LO Lagrangian, where £ = exp(2ill/f) contains the matrix of meson fields, II. It is 
useful to note that Eq. (^) can be "produced" from the LO Lagrangian of ordinary x?T by 
the substitution x — > X + P- 
The NLO Lagrangian is: 

£ 4 = L 1 <9S9S t ) 2 + L 2 <^S^S t )<^S^S t )+L 3 <(9S9S t ) 2 ) 

+L 5 (<9£d£t( x t S + + W 5 + £ f p)> 

+L 6 < X t S + Y)xf + W 6 <x f £ + £ f X> <P f £ + £V> 

+L 7 < X t S - T) x ) 2 + W 7 < X f £ - £ f X> <p f £ - £V> 

+L 8 (xWX + S f xS f x> + W 8 (p f Ex f S + £Wx) • (10) 
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The angled brackets stand for traces over the flavor indices. In the limit a — > 0, the Lj's 
above are the usual Gasser-Leutwyler (GL) coefficients of 

A word about certain log(aA) corrections is appropriate here. In the EFT formulation, 
the Lagrangian is written in terms of the most general set of operators constructed out 
of the relevant degrees of freedom that respect the symmetries of the theory. The high- 
energy physics that was integrated out enters through the unknown couplings that multiply 
these operators. Thus the low-energy constants or couplings are entirely determined by the 
high-energy scales. In x?T, the operators contain only the light meson and photon fields, 
and the low-energy constants B , f, Lj, etc., are functions of the QCD scale Aqcd- In 
particular, the Lj's are independent of the pion masses m n ~ y/mq which are associated 
with the long-distance physics. All the m q dependence of is explicitly written in the 
operators. This still holds for the m q dependence of the W%PT Lagrangian written above, 
but the same cannot be said about the lattice spacing a. It is true that an 0{a) term 
breaks the chiral symmetry in the same way as m q , and a A 2 is a soft scale associated with 
the pseudo GB mass, but 1/a is not a soft scale — it acts as the ultra-violet cutoff for 
the discrete lattice. Thus, while the low-energy constants of W%PT are expected to be 
independent of m q and acsw, they could in principal have a complicated dependence on the 
gauge coupling g (in So), which itself depends on a. However, the running of the coupling 
constant is determined in simulations by requiring that as one approaches the continuum, 
some chosen physical quantity remains fixed. Effectively, this means that the coupling g and 
the cutoff 1/a combine to give the only real scale in the theory - Aqcd ~ an d the continuum 
limit is approached smoothly. We therefore expect the Lj's and WiS to depend on Aq C d, 
and only weakly on a, the latter dependence coming from higher orders in perturbation 
theory, or involving higher powers of a which can always be expanded. 1 The parameter csw 
in the action Si will still depend on log(aAQC£>), and it might be possible to calculate these 
dependences explicitly in perturbation close to the continuum [Q, |5|. 

A related form of implicit a dependence exists in the quark masses. The quark masses 
that appear in the W%PT Lagrangian (Eqs. (Q), flTU|)) are not the same as those that appear 
in the Wilson action (Eq. ([!])). Because of the explicit breaking of chiral symmetry due to 
the Wilson term, the quark masses are not protected from additive renormalization of the 
order of the lattice cutoff 1/a. In practice one finds in simulations a "critical" line m q (a) 
on which the meson masses approximately vanish. The quark mass is then defined as the 
distance from this line: 

m q = m q - m c q (a), (11) 

and it is fh q that should be used in Eqs. (||), fllCf ). m q compensates for the large Oil/ a) 
shift in the quark masses, but it also contains positive powers of a. This is not a problem 
- redefinitions of the mass parameter of this sort only lead to changes in the W^s. The 
GL coefficients are not affected because the operators with which they are associated do 
not contain a. The re-shuffling of the W^s does mean, however, that their actual numerical 
values depend on the prescription that is used to determine m c Aa) and to define fh q . 

Note that the chiral limit cannot be taken by simply setting m q — > 0. While m c q satisfies 
M 2 (m q (a),a) = 0, there is no reason that other quantities will attain their chiral limit for 
this value of m q . This is a reflection of the fact that there really are 2 different operators 
that break the symmetry. 



We thank Paulo Bedaque and Andrew Cohen for helping us understand this issue. 
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III. APPLICATIONS 



In the following two subsections we calculate the expressions for the mass and decay 
constant of a flavor- charged meson with the flavor indices AB (A ^ B) having the same 
quantum numbers as ipBlb^A- in the calculations that follow we take x to be a diagonal 
matrix with entries (x)u = Xi-, an d use the notation xab = (xa + Xb)/2. Note that this 
notation coincides with the standard way of denoting matrix elements only for the diagonal 
entries. The same convention is used for p. It is convenient to define another matrix, 
p = x + Pi which is the combination that appears in £2- The subscript notation for p 
follows that of x an d p, except for the quantities p, n and p v which are defined below. 



A. Masses 

The mass of a flavor-charged meson in W%PT with three quark flavors is given through 
NLO by 

Mab = (^Ab)lO + (^Ab)nLO, loop + (^Ab)nLO, tree , 

(12) 

with 

{M 2 AB ) LO =Pab , (13) 
( M 1b)nlo, loop =48fvF I*ab E Rt*** lo §^ > ( 14 ) 

X=TT,r) 

24 24 8 

( M 1b)nLO, tree = ~ j2 Li ^ AB + P AB ^ * ~ J^iXABp ~ J^U (xAB + PAb) XAB 

8 24 24 

_ -pW 5 XABPAB + j^2L 6 xabX + ~p W(i (xabP + Pabx) 

24 8 

+ j^2L 8Xab + j^2W 8 XabPab , (15) 

where p n and p v are the squares of the LO masses of the two light flavor-neutral mesons, 
given implicitly by 

Pn + P v = 2p , (16) 

PttPt) = PlP2 + PlP>3 + P2P3 ■ (17) 

Here X = ^ r (x)/3 and similarly for p and p. Also, if we denote by C the flavor that is 
different from both A and B, we have 

j^AB __ PC — Pit j^AB _ PC — Prj 
pr) Ptt p-rr Prj 

In deriving the expressions for the mass in W%PT one could use a very convenient "trick" 
relating these expressions to the corresponding expressions in ordinary %PT. As mentioned 
earlier, the LO Wilson chiral Lagrangian £2, Eq. (§), can be obtained from xPT LO La- 
grangian by the simple substitution x ~ * X + P ( or X "~ > p)- Thus any quantity h(x) in 
xPT that depends only on the LO Lagrangian can be trivially reproduced in W%PT accord- 
ing to h(x) — > h(x + p). This is true for the LO and NLO loop diagrams that contribute to 
the mass. Similar results also hold for the decay constant. We provide the expressions for 
the mass in ordinary %PT in Appendix |A| for comparison. 
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B. Decay constants 



The decay constant is given through NLO by: 

fAB = (/ab)lO + (/ab)nLO, loop + (/ab)nLO, tree 

with 



(/ab)lo -/ ' 



(Jab) 



NLO, loop g 4?r 



A-BJNLO, tree 



— to lo S W + -V>b)\ log 

J i=l,2,3 J ^ 



(19) 



(20) 



Hb 



+ 



j=A,B 



+ £ ^ 



x=ir,ri 



log (fJ, A / fJ-x) l0g(/i B /^ a 



A 4 a — 



A*B — A*x 



12 4 

y (L 4 x + W^p) + - (L 5X ab + W 5 p AB ) . 



(21) 
(22) 



C. WxPT, 0(a 2 ), and improvement 

In the simplest sense, the expressions for the mass and decay constant in WxPT can be 
used to aid in taking the continuum limit. These forms provide all the linear a dependence, 
as well as non-trivial logarithms that involve a and m q . A test of these formulae would be 
to check whether they describe the a dependence better than naive extrapolations. Perhaps 
a more useful way to think about it is that with these expressions one can determine the 
GL coefficients directly from lattice data at finite a. 

What about higher orders in a? At order a 2 the picture changes qualitatively. There are 
operators in 5*2, such as ^IpD^D^tp, that do not break the chiral symmetry. This means 
that a can no longer be associated only with symmetry breaking effects, and spurion analysis 
cannot be used to constrain the a 2 operators. Nevertheless, we might still expand in e and 
S simultaneously. The LO, 0(e, 5) Lagrangian, and consequently the LO mass and decay 
constant are unchanged. At NLO, C(e 2 , eS, 5 2 ), there are several 0(a 2 ) operators that are 
added to the Lagrangian, but they are all independent of the quark masses and do not 
contain derivatives. Consequently, the only correction to the meson masses at this order is 
an additional term of the form ua 2 , where u is an unknown constant of mass- dimension 4. 
The expression for the decay constant does not receive any corrections at this order. This 
is because tree level contributions to the decay constant can only come from operators that 
contain derivatives. 

Improvement schemes (first suggested by Symanzik in Ref.0) - - using an improved 
action and improved operators — are another important tool for studying and reducing 
discretization effects. Using improved action involves adding a discretized version of the 
Pauli term in Eq. (U) to the Wilson action in Eq. (|l|) which exactly cancels the S\ term 
in the continuum action, Eq. (||]). This means that up to 0(a 2 ) the lattice theory is just 
QCD, and at low-energies we expect to be a good description. From the perspective of 
WxPT this is equivalent to saying that using an improved action sets all W{ = 0. This is of 
course not surprising: the use of improved action is meant to eliminate all 0(a) dependence 
from observables, and the Wi coefficients parameterize exactly this dependence. One should 
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not conclude from this that improving the action is enough for complete 0{a) improvement. 
As mentioned above, some dimension-5 operators that are allowed by the symmetries, and 
can therefore appear in Si, are implicitly absorbed in Sq by replacing the bare parameters of 
Sq with renormalized ones. This is a necessary step, and it is compatible with the fact that 
in improvement schemes, in addition to using an improved action, one must use improved 
operators. 

IV. PARTIALLY QUENCHED THEORIES 

W%PT is appropriate for the type of lattice simulations which are called "unquenched" . 
These are simulations in which there are 2 or 3 dynamical fermions (also called "sea quarks"), 
and expectation values are calculated of operators which are constructed from a different type 
of fermions ( "valence quarks" ) which have the same masses as the sea quarks. In most lattice 
simulations, however, the masses of the valence quarks are not taken to be the same as those 
of the sea quarks. Simulations that are done this way are called partially quenched (PQ). 
Theoretically this can be described by a QCD-like construction which includes ghosts || f7| . 
The low-energy behavior of these theories is described by PQ %PT 0, which has the same 
unknown low-energy constants as %PT for ordinary QCD ||. Thus, PQ simulations provide 
additional mass parameters that can be used to probe the theory in a larger parameter space, 
and gain better statistics in determining the GL coefficients || [I0| . It is of clear practical 
value to consider the generalization of W%PT to the PQ case. 

PQ QCD contains three different types of spin-half particles — valence quarks, sea quarks, 
and ghosts which obey Bose-Einstein statistics. There is a single ghost flavor for every 
valence quark, and they both have the same mass. The quark mass matrix for a theory with 
2 valence quarks, 3 sea quarks, and 2 ghosts is 



xPT for PQ QCD is constructed in terms of this matrix, or in terms of x which is still defined 
through Eq. (||), and the result is a Lagrangian identical to the one for ordinary x?T, but 
with an extended flavor structure and with super-traces replacing the traces. Because of 
the great formal similarity between QCD and PQ QCD, the extension of PQ %PT to PQ 
W%PT is a simple generalization of the discussion in the previous section. In particular, the 
LO and NLO Lagrangians for PQ W%PT have forms just like in Eqs. @, fllPf) , with traces 
replaced by super-traces. Further, as in the continuum PQ %PT, the low-energy constants 
L^s and WiS are exactly the same in the unquenched and PQ W^PT. It follows that one 
can use the Wilson chiral expressions for PQ theories to extract the GL coefficients. 

For completeness, we provide the expressions for the mass and decay constant for PQ 
W%PT in Appendix 0. Again, as in unquenched theories, the LO and NLO loop results in 
PQ W%PT are trivially related to the corresponding results in PQ %PT which have been 
calculated in ||. 

V. SUMMARY 

We constructed a low-energy EFT, W^PT, of the Wilson lattice action close to the con- 
tinuum. The theory extends xPT, and the perturbative framework is described in terms of 




(23) 



valence 



sea. 



ghost 
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two small parameters — the quark mass m q and the lattice spacing a. The Gasser-Leutwyler 
chiral Lagrangian (through 0(p 4 ) in xPT) was modified to incorporate all linear dependence 
on a. We applied this theory to calculate light meson masses and decay constants. The re- 
sulting expressions capture all the linear dependence on a as well as non-trivial logarithms 
that entangle a and m q . A useful application of this theory is the determination of the 
Gasser-Leutwyler coefficients of ordinary xPT from lattice simulations at small but finite a. 
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APPENDIX A: X PT RESULTS 

We present the expressions for the mass of a flavor- charged light meson in ;\RT for 
comparison with W;\PT results. As explained in the text, one can see that the LO and 
NLO loop expressions in W%PT can be obtained from the corresponding %PT results with 
the substitution \ ~~ * X + P (x "~ > /f)- Using the same notation as in Eqs. (|T3|)-(|T%|), the 
masses through NLO with three quark flavors are |§: 



with 



where 



Mab — (^Ab)lO + (^Ab)nLO, loop + (^4b)nLO, tree , (Al) 

(M 2 A b) lo = Xab , (A2) 
( M ab) nlo , 1oop = 48^5*4* E ***** log X, , (A3) 



x=ir,r] 



24 8 

( M2 Ab) NLO , tree = f2 ^ ~ ^) XABX + J 2 (2L 8 - L S ) X \b , (A4) 



+ Xv = 2 X , (A5) 
XnXv = X1X2 + X1X3 + X2X3 , (A6) 

j^AB _ XC — Xn j^AB _ XC — Xy 
Xy-Xn' V Xn-Xy' 



(Here, again, C is the flavor that is different from both A and B.) 



APPENDIX B: PQ W^PT RESULTS 

The forms of the Lagrangians in unquenched and PQ theory are the same, with an implicit 
difference in the structure of the matrices in flavor space and the replacement of traces with 
super-traces. Thus all tree contributions in both theories have the same dependence on \ 
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and p. In particular, the LO and NLO tree results are still given by Eqs. (|T3|), (|T5|) and 
Eqs. (^I]), ( p2|) for the mass and the decay constant respectively, with appropriate x an d P 
matrices for PQ simulations (the structure of p in PQ W%PT is determined by r in the PQ 
version of the Wilson action. The latter must have a structure similar to m q , eq. (|23|), that 
is needed to guarantee the exact cancellation between valence and ghost loops.). We give 
here only the NLO loop results, which are different from the unquenched expressions. 



(M 2 AB ) 



NLO, loop 



(Jab) 



NLO, loop 



1 



lQf 2 ir 2 N 



PAB ^ RxPx log(> x ) 



x=A,B ,ir,r] 



J i=l,2,3 J v 

[paD a + PbD b + (pa - Pb) 2 ] 



64tt _ 

j=A,B 

\og{p A /p B ) 



Pa — Pb 
+ ^ R xPx(pa ~ Pb) 



x=-k, t) 



\og(p A /Px) \og(p B /p x ) 



Pa — Px 



Pb — Px 



(Bl) 



(B2) 



where 



Rx 



n (pi - p* 

i=l,2,3 



II (Py - Px) ' 

y = A, B,TC,rj 



D ,. 



n (pt - Px) 

i=1.2,3 



(Pit- Px)(Pr, 



Pi 



(B3) 
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